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Abstract 

In this paper we will study the presence of extended foreground correlation features still 
present in the cleaned CMB map. It is customary to employ a galctic mask on th cleaned 
CMB map to eliminate the galactic region still contaminated even after cleaning. In this 
study we find that there are residual foreground features still present outside a galactic mask, 
which could potentially be responsible to some of the claimed large scale anomalies found 
in CMB. We defined two correlation statistics, one in pixel space and the other in multipole 
spcae. Pixel correlation statistic was used to convey a pictorial impression of the level of 
foregrounds still present in the cleaned map. The multipole correlation statistic was used 
to identify severly polluted modes in a cleaned data. Using the pixel statistic, we find that 
there are significant spatially extended features which are correlated with foregrounds. An 
interesting feature we found is a significant positive correlation of the anomalous cold spot 
region of the CMB map with synchrotron templates. When we use a multipole statistic, we 
find that the low multipoles, particularly the quadrupole is contaminated by synchrotron and 
dust. 



1 Introduction 

Cosmic microwave background radiation (CMBR) is one of the corner stones in establishing the 
inflationary Big Bang model as our current standard model of cosmology. The ACDM model or 
the cosmic concordance model based on Cosmological Principle, which assumes that the energy 
density of the universe is now dominated by dark matter and dark energy, fits the observations 
quite well [l}|3]. 

There have been many studies in the cosmology literature which found and also proposed 
various explanations to some of the large scale anomalies seen in CMB radiation. Some of these 
studies include directional and power anisotropics such as the preferred directions for universe as 
suggested by the anomalous alignment of quadrupole and octopole of CMB (4-12 and the power 



anisotropics such as the 15-35 . Since CMB is a key to understanding our universe, the presence 
of these anomalies and their origin need to be investigated thoroughly. They may lead us to new 
physics and/or provide new information about the early universe. 

An important aspect of the CMB studies is to separate non-CMB signals from the CMB, 
efficiently amidst the froth of all components registered by the detector. Even after cleaning the 



raw data, there are foreground residuals still present in the cleaned data 36-40 . The standard 



practice in CMB analysis is to mask the heavily contaminated galactic plane and other smaller 
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foreground regions and point sources from the full-sky data, by applying a foreground mask, 
before using it in estimating any cosmological parameters of interest. 

The fluctuations in CMBR, denoted as AT(fi) = T(n) — T, where T is the mean tempera- 
ture of CMB over the sky and AT(n) represent the primordial anisotropics around that mean 
temperature. We are observing CMB anisotropies as projected onto a sphere. Hence, they are 
conventionally expanded in terms of spherical harmonics, Y/ m (n), as 

oo +1 

AT(n) = Y, a i™Yi m {n) (1) 
1=0 m=-l 

where a\ m are the spectral coefficients. The spherical harmonics are suitable for expanding a 
function defined on a sphere, with n = (9, (f>) denoting the co-ordinates on the sphere, as they 
form a complete basis. The zero mean Gaussian random fluctuations of CMBR in multipole space 
have the statistical properties, (a/ m ) = and (a* lm ai> m i) = Ci6u>5 mm >, where C\ is the two point 
correlation function in multipole space, which is also called the "power spectrum and the angular 
brackets denote an ensemble (statistical) average. 

In this paper, we define two correlation statistics which can shed some light on the residual 
foreground correlations still present in the data after cleaning. The first one is defined in pixel 
space which is useful to map-out large scale foreground features / correlations hidden in the cleaned 
CMB map. The second one is defined in multipole space. The multipolar expansion coefficients, 
aim's, are widely used for variety of purposes. So this multipole correlation coefficient can reveal 
how polluted a particular multipole is. 

The CMB foregrounds are not yet well characterized. Though, thermal dust, free- free and 
synchrotron emissions from our galaxy are accepted to be the dominant foregrounds, there are 
also some anomalous non-CMB signals reported to be present in the data that contribute at 



a subdominant level. Spinning dust is one such example 15,41 



45l 57 . The presence of an 



anomalous "haze" component was found earlier in the CMB data 



46 1, though its presence was 



dismissed by WMAP team. A non-acoustic signal was found by 47] in the V — W difference 



map of WMAP's V and W bands and in 48 , an anomalous ecliptic signal was found in the 
combination map V + W — 2Q from WMAP's foreground reduced CMB maps from Q, V and 
W channels, which was speculated to be galactic straylight or zodiacal light contamination. 



Straylight contamination to CMB anisotropies were discussed in Ref. [49 , 50 . Compensating dust 
filled local voids were proposed by Ref. [5lj|52| , as an explanation to the observed CMB anomalies, 
such as the alignment of quadrupole-octopole axes, north -south ecliptic power asymmetry and 
the cold spot found in the southern galactic hemisphere of the CMB sky. In addition, the presence 
of anomalous foreground components, other than the generally recognized foregrounds viz., the 



galactic thermal dust, free- free and synchrotron emissions, are found in the WMAP data 53,54 

In view of these earlier reports, these correlation coefficients may be quite useful in assessing 
residual contaminations in the cleaned maps. 



2 Pixel and multipole space CCC's 

Correlation is one of the basic statistical measures which is employed for a wide variety of purposes. 
We try to understand the residuals in two ways, one in pixel space and the other in multipole 
space by defining suitable cross-correlation coefficients (CCC). 

Using the pixel space CCC, we try to understand the residual foreground contamination 
present in a cleaned CMB map pictorially. We do so by cross-correlating a clean CMB map 
with various foreground templates and generate an all sky correlation map. This will help us to 
understand the efficiency of the employed method for cleaning the raw data and also equips us 
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with a knowledge of the level of known and unknown contamination still present in the estimated 
CMB signal. 

The correlation coefficient in pixel space is defined as 

n CF = Y,p>e P (Cp' ~ C p ){F p , - F p ) ^ 

where C stands for the CMB signal and F stands for the foreground signal. The subscript p, is 
the pixel index. The correlation coefficient is computed in small patches across the whole sky, 
using a large pixel, p, from a low resolution HEALPi^ |55| sky discretization grid as a mask on 
both the CMB and a foreground template, which are available at higher resolution. Here p' is 
used to denote pixel indices of the high resolution CMB and foreground maps. We used CMB 
and foreground full sky maps at HEALPix resolution of -/V^e = 512 to generate the correlation 
maps at N s id e = 8. This way we have (512/8) 2 = 4096 pixels available to compute the correlation 
coefficients, robustly, in each larger pixel. 

Though, the correlation coefficient map provides a measure of the level of CMB and foreground 
correlations, one may not be able to readily infer the level of correlations due to the large pixel 
size at N a ^ = 8. In order to convey a direct impression of the level of correlations in various 
parts of the sky, we upgrade the low resolution correlation map at N S id e = 8 to N S id e = 512 and 
smooth this upgraded map with a beam FWHM of 450'. This way we will be able to identify 
correlations present in the cleaned CMB map on scales slightly larger than the pixel size of a 
HEALPix map at N si d e = 8, which is approximately equal to \JAk]YI x 8 2 x 180°/7r m 439.74 
arcmin. The smoothing does alter (decrease) the actual value of the correlation coefficient in 
each pixel, but aid us in gaining a pictorial understanding of the correlations better. 

We also define a correlation coefficient in multipole space as 

V \n c I 2 V \n F I 2 ' 
i—im im L^im I /ml 

where a/ m 's are the coefficients of spherical harmonic decomposition of either of the CMB or 
the foreground maps. We use the multipole space correlation coefficient for quantifying the 
significance of the observed level of correlation. 

These two coefficients will serve as complementary measures to understand CMB residuals. 
While the correlation maps are created to help us perceive the residual correlations visually, we 
can also make a quantitative estimate of the correlation coefficient in pixel space. However it 
is very tedious and computationally exhaustive when trying to estimate the significance of the 
correlation region by region (there are 12 x8 2 = 768 pixels in a HEALPix map at N s ide = 8). Hence, 
we use multipole CCC to make a quantitative assessment of the level of CMB and foreground 
correlations. 



3 Results 

We use IPSE cleaned map [4j[60] from WMAP seven year raw data and WMAP's seven year 
ILC map 1 61 1 as reference CMB maps. For foregrounds, we use pre-launch Planck Sky Model 
(PSM) foreground templates^ MEM seven year foreground maps from WMAP tearrj^] |l5| , the 

1 http://healpix.jpl. nasa.gov/ 

2 http: / / www.planck.fr/heading79.html 

3 http: / /lambda.gsfc. nasa.gov/ 
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total foreground (intensity) map^jfrom WMAP's five year data release produced using the pro- 
cedure described in Ref. and also the spinning dust template and the steep synchrotron like 
component obtained through an MCMC analysis by the WMAP teairj^ 56 57 . We also made 
use of the 408MHz Haslam map 58 dominated by synchrotron radiation, which is a composite 



full-sky map stitched by combining various surveys, and (2) the M8 dust model map at 94GHz 



of Ref. 59 



All the maps considered here are smoothed to a common 1° resolution. Some of the maps used 
here were upgraded to HEALPix resolution 9, if they are available at other resolutions. WMAP 
MEM templates are available at iV^e = 128 and WMAP MCMC foreground maps are available 
at N side = 64. 

3.1 Pixel correlation studies 

We first present our results from using pixel space CCC. As an initial step towards these cor- 
relation studies, we first correlate the cleaned CMB maps obtained using the IPSE and ILC 
procedures. The CMB maps and their correlation map are shown in Fig. [I]. As is expected, 
we find that these two maps are highly correlated, as can be easily inferred visually and the 
accompanying scale below the figures. Except for some regions in and close to the galactic plane, 
the CCC is close to one. 

We next discuss the spatial correlations between CMB and foreground maps. Here we cross- 
correlate full-sky maps. While extracting the CMB signal, one generally removes the most con- 
taminated region. The KQ85 galactic mask, used in CMB analysis of WMAP's seven year tem- 
perature data, is shown in Fig. |2j. We expect significant correlation between the foregrounds and 
CMB signal in this region. However, outside the masked region, we expect that the correlation 
would be negligible. 

In Fig. [3] and Fig. [4j, the correlations between IPSE cleaned CMB map [4,60 , obtained 



from WMAP seven raw data, and foreground templates of galactic thermal dust, free- free and 
synchrotron emissions, from pre-launch PSM and WMAP seven year MEM foreground templates, 
are shown. In Fig. [5], correlations with spinning dust and steep synchrotron like component 
from an MCMC analysis by the WMAP team are presented. 



4 https: / /www. cfa.harvard.edu/~adeolive/gsm/index. html 
5 http: / /lambda.gsfc. nasa.gov/ 
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Figure 2: The KQ85yr7 galactic mask used by WMAP team in their CMB temperature analysis 
to eliminate foreground contaminated regions in the cleaned data. 



In Fig. [6], the correlation between the IPSE cleaned map with Haslam's 408 MHz map, 



dominated by synchrotron emission 58 and the M8 dust template from Ref. 59 are shown. Both 
these maps, and extinction corrected Ha map as free- free emission template, are used as priors in 
generating the MEM foreground templates at WMAP's frequencies. Finally, the correlation maps 



of IPSE cleaned CMB map, and the total intensity foreground maps from Ref. 11 produced from 
WMAP's five year data are shown in Fig. [7j. 

In the figures JsJ , jjjj, [l0| , 11 and [l2| , the same correlation for the seven year ILC ma F0|6l| 
from WMAP team are shown in the same order. 

We have also extended our analysis to a CMB map generated using the recently introduced 
kurtosis minimization, instead of variance minimization in pixel space [62] . The clean map ob- 
tained so was referred to as Gauss map. The spatial correlations between the Gauss map and 



the pre-launch PSM and seven year WMAP's MEM foreground templates are shown in Fig. 13 



and Fig. 14 , respectively. 

One could readily notice extended spatial correlations across the CMB sky. Some of these 
features appear very prominent. Some of our immediate observations are 

• Correlation with PSM dust template near the northern galactic poles. 

• Large positive correlation of the cleaned CMB maps with synchrotron maps roughly in the 
direction of the southern cold spot. The cold spot from the IPSE cleaned CMB map from 



WMAP seven year data is shown in Fig. 15 



Similar spatial correlation in the direction of the cold spot (roughly) are seen with spinning 
dust and steep synchrotron like component templates from WMAP's MCMC analysis. 

Significant correlations with 408 MHz Haslam map, especially the spur like region near the 
south pole which could be related to the redder regions along the ecliptic found by Ref. [48] 
using the V + W — 2Q combination map from foreground reduced WMAP CMB band maps. 



J http: / /lambda.gsfc. nasa.gov/ 
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There are also correlations with dust M8 model template from Ref. [59] around the north 
pole. 



The ringer like structure in the correlation maps with the total foreground (intensity) maps 
[111 on the left side of the images. These finger like extensions, of which the southern cold 



spot is also a part of, were pointed out by WMAP team in their seven year analysis 63 
The positive correlation in the direction of cold spot fades away in the correlation map with 
total foreground map at W— band. So, it suggests that it is originating from a foreground 
component(s) at low frequencies, which could be synchrotron. 

• Interestingly enough, we don't find such an identifiable correlation for the large cold spot, 
present slightly to the right from the center of the CMB map in galactic co-ordinates, with 
any of the foregrounds considered here. 

There may be many features which we may not have identified here. This spatial correlation 
analysis can help us understand the residual foregrounds better. This correlation analysis using 
temperature data can also be trivially extended to include correlations with Q and U polarizations. 

We emphasize that these are mere observations and an eventual speculation with known 
anomalous features, reported in CMB so far. A consistent further analysis needs to be done 
before making any strong claims. 



The cold spot [22-24 64 -75 correlation feature is an interesting finding. We do a preliminary 
analysis of this feature here. We generate a sky mask of 10° x 14° size centered in the direction of 
the cold spot. This cold spot is roughly located at (0, (p) = (147°, 209°) [22]. We apply this patch 
simultaneously on the synchrotron template and the CMB map, and compute the correlation 
from that region. We find that 1Z P = 0.4368. In order to estimate it's significance, we applied the 
same mask on 800 simulated cleaned IPSE maps and corresponding pure CMB maps, both at one 
degree smoothing, separately, and correlated with the synchrotron template from PSM. We also 
did a similar correlation analysis in the direction of cold spot using 600 simulated ILC cleaned 
maps and the corresponding pure maps. The results for both ILC and IPSE cleaned simulated 
CMB maps and pure maps as correlated with synchrotron, free-free and thermal dust emission 



templates, respectively, from PSM are shown in Fig. 16 , Fig. |17| and Fig. |18 . The normalized 



frequency counts of the CCC from simulated pure CMB data are peaked at zero, as expected. 
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Figure 7: The IPSE cleaned CMB map is correlated with total intensity foreground maps by 
using the procedure of Ref. from five year WMAP data. From top to bottom rows are the 
total foreground (intensity) maps at WMAP's K, Ka, Q, V and W bands. Notice the finger like 
extensions on the left side of the images, extending from the galactic plane towards the southern 
hemisphere. The cold spot is pointed out to be a part of this finger like structure by the WMAP 



team in their "CMB anomalies" paper 63 usirjg? seven year data. 
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Figure 15: The cold spot of the CMB map from the IPSE cleaned map from using WMAP seven 
year data. 



The CCC histogram plot of the frequency counts for the correlation of a foreground and 



simulated CMB maps (pure or clean) in the direction of the cold spot are shown in Fig. 16 
We find that the histogram plots generally broaden for the cleaned maps. The pixel correlation 
coefficient, 1Z P , for the data are shown for comparison as coloured dots in those plots. In ILC 
map, the cold spot correlation histograms for pure maps and clean maps are similar to each other. 
From this preliminary analysis, the cold spot correlation with dust and synchrotron are outside 
2a CL and the correlation with free- free emission is insignificant. 

3.2 Multipole space correlations 

In this section we quantify the level of foreground residual contamination using multipole space 
CCC. We cross correlate the IPSE cleaned seven year map, and WMAP's seven year ILC map, 
with the three principal foreground maps from PSM in multipole space following Eq. [3j. We 
also compared the data correlation coefficients using 800 simulated IPSE cleaned maps and 600 
simulated ILC cleaned maps, and also using the corresponding pure maps. The results are 



shown in Fig. 19 , Fig. [20] and 21 , for multipole space correlation with dust, free-free and 
synchrotron emissions, respectively. A la confidence level with the mean correlation estimated 
from simulations is also plotted. We could readily find that some of the normalized correlation 
statistic values are outside the la CL. The farther these data CCC's are from the mean, the more 
polluted they are. Of course, these are correlations with foregrounds and full-sky cleaned CMB 
maps. As mentioned earlier, a galactic mask is used to avoid foreground induced errors in CMB 
analysis. How these statistics change with masking is not known at present. But some studies 
use the spectral coefficients(a/ m 's) from full-sky maps, such as the directional anisotropy studies. 
Hence, there this information from multipole CCC may be useful. 
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Figure 16: The correlation of the cold spot w ^h PSM synchrotron map at 23GHz with IPSE 
(top) and ILC (bottom) seven year cleaned maps. 
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Figure 17: The correlation of the cold spot reg^n with PSM free- free template at 23GHz with 
IPSE7yr (top) and ILC7yr (bottom) cleaned maps. 
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CCC of cold spot and PSM Dust at 94GHz (800 IPSE cleaned maps) 

CCC from IPSE7yr map • 



0.25 



c 

50 



0.2 



0.15 



0.1 



0.05 



-1 



-0.5 



0.5 




Figure 18: Correlation statistics for the cold S P,°^ region with dust template at 94GHz of PSM 
with IPSE (top) and WILC (bottom) seven year cleaned maps. 
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Figure 19: Multipole correlation of IPSE (top) and ILC (bottom) cleaned maps with PSM dust 
map at 94GHz.The extend of the vertical bars denote the la CL as estimated from simulations. 
This graph reveals us the most polluted multipoles in th CMB map. Note that these correspond 
to full-sky CMB and foreground correlations 
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Figure 20: Same as Fig. 19 , but correlated with PSM free-free emission map at 23GHz. 
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Figure 21: Multipole correlation statistics from correlating IPSE (top) and ILC (bottom) cleaned 
maps with 23GHz synchrotron map of pre-launch Planck Sky Model. 
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4 Conclusions 



In this paper, we defined two correlation statistics, one in the pixel space and the other in the 
multipole space to understand the foreground residuals still present in a cleaned CMB map. 
We were particularly interested in determining the presence of such residuals outside the KQ85 
mask used in CMB temperature analysis. Discarding a portion of the CMB sky near heavily 
contaminated galactic plane by applying a mask is a preventive measure in CMB analysis, to avoid 
foreground induced errors. We find many extended spatial correlation features across the CMB 
sky with various foreground templates that are assumed to contaminate the CMB measurements. 
We studied the correlations with both the three principal foregrounds viz., synchrotron, dust and 
free-free galactic emissions, and the subdominant foregrounds such as the spinning dust and the 
steep synchrotron like spectral index. For these we used the foreground maps which are available 
in the public domain. 

One interesting feature we identified is that there is a large positive correlation between 
the CMB map and the synchrotron template roughly in the direction of the cold spot situated 
near the galactic south pole. A preliminary investigation of this observation was also done by 
correlating that part of the sky by applying a 10° x 14° mask in that direction on both the pure 
maps and cleaned maps with synchrotron template. We also did similar correlation analysis with 
dust and free-free foreground templates. The multipole space correlation revealed how the various 
multipoles are contaminated by foregrounds. We find that many of the multipoles estimated from 
full sky cleaned map are contaminated, as shown by simulations with mock cleaned maps. The 
quadrupole is found to be significantly contaminated by the synchrotron component. However, 
since the multipole CCCs are computed by correlating full sky maps, directly inferring the effect 
of the residual foreground features found here may not be obvious. This resutls of this study to 
some extect point outs the ambiguity in te foreground separation and possible contribution of 
residuals to any of the anomlaies found in CMB data. 
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References 

[1] Riess A. G. et al., 1998, The Astronomical Journal, 116, 1009 
[2] Perlmutter S., et al, 1999, ApJ, 517, 565 
[3] Larson D., et al, 2011, ApJS, 192, 16 

[4] Tegmark M., de Oliveira-Costa A. and Hamilton A. J., Physical Review D,68, 123523 

[5] de Oliveira-Costa A., Tegmark M., Zaldarriaga M. and Hamilton, A., 2004, Physical Review 
D, 69, 063516 

[6] Copi, C. J., and Huterer D. and Starkman G. D., 2004, Physical Review D, 043515 
[7] Ralston, J. P. and Jain P., 2004, 13, 1857 
[8] Souradeep, T. and Hajian, 2004, A, Pramana,62,793 
[9] Hajian, A. and Souradeep, T. and Cornish, N., 2005, ApJL,618,L63 
[10] Land, K. and Magueijo, J., 2005, Physical Review Letters, 95, 071301 




27 



[11] de Oliveira-Costa, A. and Tegmark, M., 2006,Physical Review D, 74,023005 

[12] Copi, C. J. and Huterer, D. and Schwarz, D. J. and Starkman, G. D., 2006,MNRAS,367,79 

[13] Samal, P. K. and Saha, R. and Jain, P. and Ralston, J. P., 2008, MNRAS, 385, 1718 

[14] Samal, P. K. and Saha, R. and Jain, P. and Ralston, J. P., 2009,MNRAS,396,511 

[15] Bennett C.L.., et al., 2003,ApJS, 148, 1 

[16] Efstathiou, G., 2003,MNRAS,346,L26 

[17] Efstathiou, G., 2004,MNRAS, 348,885 

[18] Copi, C. J., Huterer, D., Schwarz, D. J. and Starkman, G. D.., 2007, Physical Review 
D,75,023507 

[19] Eriksen, H. K. and Hansen, F. K. and Banday, A. J. and Gorski, K. M. and Lilje, P. B., 
2004, ApJ, 605, 14 

[20] Eriksen, H. K. and Banday, A. J. and Gorski, K. M. and Hansen, F. K. and Lilje, P. B., 
2007, ApJL, 660, L81 

[21] Hoftuft, J. and Eriksen, H. K. and Banday, A. J. and Gorski, K. M. and Hansen, F. K. and 
Lilje, P. B., 2009,ApJ,699,985 

[22] Vielva, P. and Martinez-Gonzalez, E. and Barreiro, R. B. and Sanz, J. L. and Cayon, L., 
2004, ApJ, 609, 22 

[23] Cruz, M. and Martinez-Gonzalez, E. and Vielva, P. and Cayon, L., 2005,MNRAS, 356,29 

[24] Cruz, M. and Tucci, M. and Martinez-Gonzalez, E. and Vielva, P., 2006, MNRAS, 369, 57 

[25] Ackerman, L. and Carroll, S. M. and Wise, M. B., 2007, Physical Review D, 75, 083502 

[26] Groeneboom, N. E. and Eriksen, H. K., 2009, ApJ, 690, 1807 

[27] Hanson, D. and Lewis, A., 2009, Physical Review D,80,063004 

[28] Hanson, D. and Lewis, A. and Challinor, A., 2010, Physical Review D, 81, 103003 

[29] Land, K. and Magueijo, J., 2005,Physical Review D, 72, 101302 

[30] Gurzadyan, V. G. and Starobinsky, A. A. and Kashin, A. L. and Khachatryan, H. G. and 
Yegorian, G.,2007, Modern Physics Letters A, 22, 2955 

[31] Kim, J. and Naselsky, P., 2010, ApJL, 714, L265 

[32] Kim, J. and Naselsky, P., 2010, Physical Review D, 82,063002 

[33] Gruppuso, A. and Finelli, F. and Natoli, P. and Paci, F. and Cabella, P. and de Rosa, A. 
and Mandolesi, N., 2011, MNRAS, 411, 1445 

[34] Aluri, P. K. and Jain, P., 2012, MNRAS,419,3378 



[35] Ben-David, A. and Kovetz, E. D. and Itzhaki, N.,2011, arXiv: 1108. 1702 



[36] Slosar, A. and Seljak, U., 2004, Physical Review D, 70, 083002 

28 



[37] Abramo, L. R. and Sodre, Jr., L. and Wuensche, C. A., 2006, Physical Review D, 74, 083515 
[38] Rakic, A. and Rasanen, S. and Schwarz, D. J., 2006,MNRAS,369, L27 
[39] Gruppuso, A. and Burigana, C, 2009, JCAP, 8, 4 

[40] Aluri, P. K. and Samal, P. K. and Jain, P. and Ralston, J. P., 2011, MNRAS,414,1032 

[41] Kogut, A. , Bulletin of the American Astronomical Society,1996, 28, 1295 

[42] de Oliveira-Costa, A. and Kogut, A. and Devlin, M. J. and Netterfield, C. B. and Page, L. A. 
and Wollack, E. J., 1997, ApJL, L17 

[43] Leitch, E. M. and Readhead, A. C. S. and Pearson, T. J. and Myers, S. T., 1997, ApJL,486, 
L23 

[44] de Oliveira-Costa, eta 1., 2002, ApJ, 567, 363 

[45] Dobler, G. and Finkbeiner, D. P., 2008,ApJ,680,1235 

[46] Dobler, G. and Finkbeiner, D. P., 2008, ApJ, 680,1222 

[47] Jiang, B.-Z. and Lieu, R. and Zhang, S.-N. and Wakker, B., 2010,ApJ, 708,375 

[48] Diego, J. M. and Cruz, M. and Gonzalez-Nuevo, J. and Maris, M. and Ascasibar, Y. and 
Burigana, C, 2010, MNRAS, 402, 1213 

[49] Burigana, C. and Maino, D. and Gorski, K. M. and Mandolesi, N. and Bersanelli, M. and 
Villa, F. and Valenziano, L. and Wandelt, B. D. and Maltoni, M. and Hivon, E., 2001, A&A, 
373, 345 

[50] Burigana, C. and Gruppuso, A. and Finelli, F., 2006, MNRAS,371, 1570 

[51] Inoue, K. T. and Silk, J., 2006, ApJ, 23 

[52] Inoue, K. T. and Silk, J., 2007, ApJ, 664, 650 

[53] La Porta, L. and Burigana, C. and Reich, W. and Reich, P., 2008, A&A, 479, 641 

[54] Bonaldi, A. and Ricciardi, S. and Leach, S. and Stivoli, F. and Baccigalupi, C. and de Zotti, 
G., 2007, MNRAS, 382, 1791 

[55] Gorski K. M., Hivon E., Banday A. J., Wandelt B. D., Hansen F. K., Reinecke M., and 
Bartelmann M., 2005, ApJ, 622, 759 

[56] Gold, B. eta 1., 2009,ApJS,180, 265 

[57] Gold, B. et al., 2011,ApJS, 192, 15 

[58] Haslam, C. G. T. and Salter, C. J. and Stoffel, H. and Wilson, W. E., 1982, A&AS, 47, 1 
[59] Finkbeiner, D. P., 2003, ApJS, 146,407 

[60] Saha, R. and Prunet, S. and Jain, P. and Souradeep, T., 2008, Physical Review D, 78, 023003 
[61] Hinshaw G., et al, 2007, ApJS, 170,288 
[62] Saha, R., 2011,ApJL, 739,L56, 

29 



[63] Bennett. C. L., et al., 2011, ApJS, 192, 17 

[64] Rudnick, L. and Brown, S. and Williams, L. R., 2007, ApJ, 671, 40 

[65] Cembranos, J. A. R. and de la Cruz-Dombriz, A. and Dobado, A. and Maroto, A. L., 
2008,JCAP, 10, 39 

[66] Cruz, M. and Martinez-Gonzalez, E. and Vielva, P. and Diego, J. M. and Hobson, M. and 
Turok, N., 2008, MNRAS, 390,913 

[67] Bridges, M. and McEwen, J. D. and Cruz, M. and Hobson, M. P. and Lasenby, A. N. and 
Vielva, P. and Martinez-Gonzalez, E., 2008,MNRAS, 390,1372 

[68] Masina, I. and Notari, A., 2009,JCAP, 2, 19 

[69] Bernui, A., 2009,Physical Review D, 80,123010 

[70] Smith, K. M. and Huterer, D., 2010, MNRAS, 403, 2 

[71] Naselsky, P. D. and Christensen, P. R. and Coles, P. and Verkhodanov, O. V. and Novikov, 
D. I. and Kim, J., 2010, Astrophysical Bulletin, 65, 101 

[72] Solov'ev, D. I. and Verkhodanov, O. V., 2010,Astrophysical Bulletin, 65,121 

[73] Bremer, M. N. and Silk, J. and Davies, L. J. M. and Lehnert, M. D., 2010, MNRAS,404, 
L69 

[74] Granett, B. R. and Szapudi, I. and Neyrinck, M. C, 2010, ApJ, 714,825 

[75] Vielva, P. and Marti and Nez-Gonzalez, E. and Cruz, M. and Barreiro, R. B. and Tucci, M., 
2011, MNRAS,410, 33 

[76] Pavan K. Aluri, "Large scale anomalies in the Cosmic Microwave Background Radiation", 
Thesis submitted to Department of Physics, IIT Kanpur, December, 2011 



30 



